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Momentum-energy tensor associated to the quasiparticles in anisotropic 

superconductors 
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From a Lagrangian density for the Bogoliubov de Gennes equations in anisotropic superconduc¬ 
tors, we find the momentum-energy tensor associated to the quasiparticles of the system. For this, 
we make infinitesimal translations on both space and time and we use the Noether’s theorem. We 
prove that beyond to the usual terms associated to the electron-hole dynamic in electromagnetic 
potentials, appears terms that involves the pair potential and that are obtained from the coupling 
of the electron-like and hole-like quasiparticles. 

PACS numbers: 


I. INTRODUCTION 

The Bogoliubov-de Gennes (BdG) equations describes 
the bahaviour of the elementary excitation of a super¬ 
conductor system fli- Since there is a Lagrangian 
density that allows one to derivate these equations and 
through of the Noether’s theorem for infinitesimal trans¬ 
formations leaving the action invariant Q, it is possible 
to find the conservation’s laws associated to the super¬ 
conductor system. In this paper we build the momentum 
- energy tensor from infinitesimal transformations both 
in space and time. 


and the external potencial, ^ is the chemical potencial 
and A(r,r') is the pair potencial associated to quasipar¬ 
ticles and it is coupled to the electron and hole compo¬ 
nents. From the lagrangian density, 

if} 

^ - '0+(r,f)V'(r,t)] 

- t)]+a^P{r, t)ip{v, t)) + 

- j (3) 


MOMENTUM-ENERGY TENSOR FOR 
ANISOTROPIC SUPERCONDUCTORS. 

The Bogoliubov - de Gennes equations describe the 
quasiparticles in a superconductor are given by: 

=H'ipir,t) = 

Ho{r,t)il;{r,t)+ J dr (1) 

With, 


The Bogoliubov - de Gennes (BdG) equations are de¬ 
rived by means of the principle of least action Q . The 
Noether theorem allows us to find both a set of equa¬ 
tions and conservation laws from transformations that 
leave the action invariant. If an infinitesimal transfor¬ 
mation of p parameters that leaves action unchangeable, 
then there are p quantities conserved Q . If the motion 
equation for the fields is obtained from the princi¬ 
ple of least action and moreover the transformations in 
space, time and the fields leaves the action invariant. The 
Noether theorem yields: 


V'(r,t) = 


i(r,t) 


v{r,t) 

U(r,t) =(C7(r,t) - 7f(r,t) = 


A(rr') = f 

I A*(r,r') 0 


^ TT^ (t '' 
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7re(r,t) 0 
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dt 


V-J = 0, 


( 4 ) 


p = L5t — 


dC 




( 5 ) 


where 'ip{r,t) is the wavefunction of the quasiparti¬ 
cle with two componets electron u{r,t) and hole v{r,t), 
whereas 7fe(r, t) = {—ihV — eA) and rthix, t) = {—ihV + 
eA) depletes the electron and hole momentum respec¬ 
tively; A(r, t) is the magnetic vector potential, U (r, t) is 
the scalar potencial that includes both the Hartree-Fock 


Ji = C5xi - ^ ^ - S^kj (6) 

Where p and J correspond to the charge and current 
densities associated to the infinitesimal transformation. 
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Moreover when Eq. [6] is integrated over all space, one SPATIAL TRANSLATION-CONSERVATION OF 
obtciinsi LINEAR, IVIOIVIENTUIVI 



dQ{r,t) 

dt 


(7) 


Now we can do an infinitesimal translation in coordi¬ 
nates 


Where Q is a generalized charge and is itself the con¬ 
served quantity. 


x^=Xi + 5xi i = 1,2,3. (11) 

By using the Eq. [5] for each tspatial component x, y, 
and z we get a charge density associated given by: 


TEMPORAL TRANSLATION-CONSERVATION 
OF THE ENERGY 


Let us do an infinitesimal transformation in time, in 
this way: 


t =t + 5t. ( 8 ) 

The variations respect to the fields and spatial coor¬ 
dinates are vanished.By plugging the lagrangian density 
Eq. [3] into Eq. [®], we obtain: 

Po(r,i) =^[7r(r)'!/'(r,t)]+cr^ • 7r(r)V^(r, t) 

+ii'^{Y,t)V{T)'\l){Y,t) + / dr'ip~^{r,t)A{r,r')ip(r,t) 

(9) 

This expression can be interpreted as an energy density 
that is shown as the sum of the energy density associated 
with the electron and hole components plus the energy 
density that couples electrons and holesdepending of the 
pairs potential A. 

By using Eq. [7] and recalling that this transformation 
leaves the action invariant, we have J podV = ^ (E) = 
0, which implies that the expectation value of energy is 
a constant of motion. Now let us examine what happens 
with the density of current in this case. From Eq. [5] 
for the current density and again using the lagrangian 
density Eq. [0, we have 


— 


= i \[—ihazdj — Hijj 

+ - eAji)'4’ j 


1,2,3. (10) 


Where = v is the velocity of group of the particle. 
From this equation is observed that J is an energy flux, 
because it can be written as the quasiparticle density of 
charge times its group velocity. 


Pi = iP+TTitjj = - eA^a^)il; ( 12 ) 

In the same way. Using Eq. [6] the current density 
associated to pi is given by: 


Ji,o “2 L 




= -[[(-f?i(T^5j - eAji)^j;]^{-ihdii - eAia;,)'>p 
+ [{-iTidil - eAiaz)ii^\^{-i‘hazdj - eAjl)'ii)\ 


(13) 

(14) 


We obtain that Eq. m corresponds to the density of 
linear momentum, whereas Eq. [ [35] corresponds to the 
density of current of momentum, that is proportional to 
the density of momentum and to the velocity of group of 
the quasiparticle. If the action is invariant under these 
translations, one can integrate over all space the density 
of momentum and it yields 

S / = 4 {,.) (15) 

Therefore the momentum associated to the quasipar¬ 
ticles is conserved. The density of momentum can be 
written as pi = Pei + Phi, where Pet and phi are the mo¬ 
mentum density associated to electron and hole compo¬ 
nents, respectively. So we have 


Pei = u*{-ihdi - eAi)u ] phi = v*{-ihdi + eAi)v. 

(16) 

In the case that action is not invariant under trans¬ 
formation Eq. |12| , there will appear sources proportional 
to the scalar potential gradient, magnetic vector poten- 
cial and the pairs potential. All the equations we have 
derived so far from the space-time transformations are 
summarized in a tensorial form : 


/ Po 

Joi 

J 02 

Jo3 \ 

Pi 

Jll 

J 12 

Jl3 

P2 

J 21 

J 22 

J 23 

\ P3 

J 31 

J 32 

J 33 J 


Where the term Too corresponds to the density of en¬ 
ergy of the quasiparticle, whose explicit form is given in 
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Eq. [9], the terms of the form Toi i = 1, 2, 3 are related 
with the the density of current of energy in each com¬ 
ponent that have been calculated in Eq. m- Now the 
terms of T^q represents the density of linear momentum 
in the direction i given by Eq. |12| and finally the terms 
of the form Ty corresponds to the component j of the 
momentum current density is given by Eq. [35] . In con¬ 
clusion, the tensorial form for the energy and momentum 
equations written as: 


^ (Pi) = I (diA^Vj + (22) 

Since v = clearly Eq. [22] does not match with 

the Lorentz force definition. Here repeated indexes sum 
is considered and i = 1,2,3 are the coordinates. Let us 
find the temporal evolution of the potential vectorial A : 


-^=0 u,z; = 0,1,2,3. (17) 

dxu 

When a spatial or temporal infinitesimal transforma¬ 
tion does not leave the lagrangian density invariant an 
associated source appears. In classical mechanics, when 
a potential does not depend on the position the total 
force is zero and is a conserved quantity. In contrast, if 
the potential does depend of the position the quantity is 
not conserved and the source is the total force. In the 
next section we will see that the langrangian Eq. [3] is 
not invariant if for instance, it contains a potential vector 
dependence £(A(a;)). 


EXAMPLE: LORENTZ FORCE AS A SOURCE 
FROM A LAGRANGIAN DEPENDING ON THE 
POTENTIAL VECTOR 


A particle with mass m and electric charge q = —e in 
presence of a magnetic field can be described by: 

= (18) 

where p = —iKV. A magnetic field is generated by 
a potential vector A. From Eq. (TB] one obtains the 
continuity equation Eq. |3] with. 


J = 


2m - 


tjj~^ — qA^ 'll; — Ip -f qA^ 7p~' 


(19) 

The temporal evolution of an operator O is described 


by 



( 20 ) 


By definition the conjugate momentum is written as, 
TT = p — eA and the velocity is given by: 


X = V = 


p — eA 

m 


TT 

m 


p — eA 


( 21 ) 


therefore Tr/m plays the role of velocity and the tem¬ 
poral evolution of tt will give rise to the quantum Lorentz 
force, rather than the temporal evolution of p. In fact. 



(23) 



(24) 



{Po + • {pi + (lA 

2m 



(25) 


we can use the fact that [a, = [a, b]b + b[a, 6]. There¬ 

fore, straightforwardly 


(^7 = ^ 


A\(pj+qAj'j (pj+qAj'j 


+ (pj + A’-, ^Pj + qAjj 


(26) 


Again ones associates (jjj + qA^ = mvj = TTj. There¬ 


fore 


(iA = i 

dt\ / 2 


A 


, (^pj+qAj 


Vj + Vj 


A\ (pj + qAj 

(27) 


Evidently [A®, A-l] = 0. One ends up with: 


dt\/ 2 


A* 


Pi 


V-j + V-j 


A* 


Pj 


(28) 


Easily one shows that [A®,pj] = ihdjA^. Then, 


^ y + VjdjA^'^ (29) 

Results from either Eq. [21] and Eq. [22] give us a 
direct way to calculate the temporal evolution of the con¬ 
jugate momentum. It means. 



Plugging results from Eq. [25] and Eq. [22] yields. 
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I (r) = I ( (da^ - d,A^) %+ V, - dd^) ) 

(31) 

rewritting in terms of the magnetic field, 

= (32) 

That yields directly the definition of the Lorentz force: 

|w = |(vxB-Bxv) (33) 

it has been showed that the Lorentz force is derived 
from the canonical momentum tt rather than its conju¬ 
gate one p. 

The current equation Eq. m can be written is terms 
of velocity as: 


Constant vectorial potential 

In this case, one has that the langrangian is invariant 
under spacial translation and the associated charged is 
conserved. In other words 

= 0 (40) 

where the energy-momentum tensor is written explic- 
ity as 

(41) 

Recalling for the momentum-component: 

r° = di-ip — ihdid'4’ (42) 


J = 


1 

2m - 


— qA^ Ip — Ip -f d 


(34) 


J = 


2m - 


ip'^vip — {vtpd ^ 


(35) 


and from Eq. 


ft (Pi) = 2 + Vjdd^) = 

f dri/)+ (^iA^Vj + VjdiA^'j tp (36) 

= J dr [tp'^Vjtp + {vj'ip~^) tp] diA^ = f drPdiA^ 


Analogously, the lagrangian density can be written as: 


Position-dependence of the vector potential. 

by performing an spacial translation along the i dirrec- 
tion, it yields that 


J d^x dd + ~ = ~ J d^xduTp 


Sxi (43) 


The vector potential is non invariant under spatial 
translation. It gives: 


[C{x^ + SxA - Cd)] 


^dd^Sx, 

5Ai 


(44) 


From Eq. [55] one knows that SCj5A^ = gj*, therefore 


= Y - didd ~ ^ (37) 

and 

-Co = Y [dditp - didd\ - dAd dd (38) 

where the hamiltonian density is defined as 

^ sjc 

H = (ntp) (ntp) and Cq is the free lagrangian den¬ 
sity (lagrangian density with A = 0.) 

The lagrangian density previously written is split in 
this convinient way: 


C = Co + 3-A (39) 

and J is given by Eq. [35]. Let us considera couple of 
cases of interest: whether or not the vectorial potential 
has a dependence with position. 


J d^x [Cd + SxA - Cd)] = - J d^xduTp 


f dP^xqJ^did = -f d^xduTp 


= - J d^xdoT°Sxi = 


In conclusion. 


d = i(p- eA.) = A (p) _ e4r (A 


dt dt 


dt 


4 (p) = |(vxB-Bxv)- 4 (A 


And finally: 


Y (P*) = doreut. - I {dd^V^ + 


(45) 

(46) 

(47) 

(48) 

(49) 
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CONCLUSIONS 


From infinitesimal transformations both in space and 
time we have found the Energy-Momentum tensor associ¬ 
ated to quasiparticles for an anisotropic superconductor. 
In the case where the transformations leave the action 
invariant , we find that the energy and linear momentum 
of quasiparticles are the quantities conserved. 
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